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1. Introduction 

In the Standard Model, the flavor-changing neutral current process b — XJ^l^ only occurs at 
the one-loop level and is therefore sensitive to new physics. In the kinematical region where the 
lepton invariant mass squared is far away from the cc-resonances, the dilepton invariant mass 
spectrum and the forward-backward asymmetry can be precisely predicted using large nih expan- 
sion, where the leading term is given by the partonic matrix element of the effective Hamiltonian 
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-^eff = -^V*V,tl^Q{^)Oi{^). 
v2 ,-=1 



(1.1) 



We neglect the CKM combination V*,.Vub and the operator basis is defined as in [jl|]. In [|2|] we 
published the first analytic NNLL calculation of the high region of the matrix elements of the 
operators 

Ol = {sLy^T"cL){cLfT%L), 02 = {sLy^CL){cLy^bL) , (1.2) 

which dominate the NNLL amplitude numerically. Earlier these results were only available analyt- 
ically in the region of low q^ 01- Using equations of motion the NNLL matrix elements of the 
effective operators take the form 



{se+r\0i\b)2. 



loops 



_^\2 
471/ 



tree +Fr{09) tree 



(1.3) 



where Oj = e/ g',mh{sLO^''bR)F^y and O9 = / gi{sanbL)Wy^^) 



2. Calculations 




Figure 1: Diagrams that have to be taken into account at order a^. The circle-crosses denote the possible 
locations where the virtual photon is emitted (see text). 



The diagrams contributing at order are shown in Figure We set m., = and define 



s = ^ and z = ^ , 



(2.1) 
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where q is the momentum of the virtual photon. After reducing occurring tensor-like Feynman 
integrals [§] the remaining scalar integrals can be further reduced to master integrals using integra- 
tion by parts (IBP) identities [||]. Considering the region s > 4z, we expanded the master integrals 
in z and kept the full analytic dependence in s. 

For power expanding Feynman integrals we use a combination of method of regions and 
differential equation techniques [||, Consider a set of Feynman integrals Ii,...,!^ depending on 
the expansion parameter z and related by a system of differential equations obtained by differenti- 
ating la with respect to z and applying IBP identities: 

d \ - 

-r^a = l^Kjilii +ga, (2.2) 
az p 



where ga contains simpler integrals which pose no serious problems. Expanding both sides of ( |2.2| ) 
in £, z and Inz 

la = I/^f)£V(lnz)^ h„p='£h%e'z\ ga = E^iif £V(ln^)^ (2.3) 

'J,k i-J ij,k 

and inserting ( |2.3[ ) into ( ^^ we obtain algebraic equations for the coefficients I'^J'P 

p i' r 

This enables us to recursively calculate higher powers of z once the leading powers are known. In 



practice this means that we need the I^^f^ and sometimes also the /^'f as initial condition to (E^). 
These initial conditions can be computed using method of regions. A non trivial check is provided 
by the fact that the leading terms containing logarithms of z can be calculated by both method of 



regions and the recurrence relation ( [2.4| ). 



The summation index j in ( [2.3| ) can take integer or half-integer values, depending on the 
specific set of integrals /«. In order to determine the possible powers of z and ln(z) we used the 
algorithm described in [^]. A given D-dimensional L-loop Feynman integral I{z) reads in Feynman 
parameterization 

/ • \i . N TjN-(L+\)D/2 



where U, Fi and F2 are polynomials in jc„. Using Mellin-Bames representation (2.5) can be cast 
into the following form 

/(.) ^(-if {^y^jy^^n-sns+N-LD/2) 

X fd'^xSil-l^ x„)U''-<^'^+'W2f^F^-^-n+ld/2_ ^2.6) 

n=\ 

By closing the integration contour over s to the right hand side the poles on the positive real axis 
turn into powers of z- If we apply the technique of sector decomposition [ [l0| ] to ( |2.6| ) we end up 
with terms of the following form 

^Ik ^lJk^2Jk ' (2.7) 
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where Uik, Fij^ and F2jk contain terms that are constant in 7. From ( p77[ ) we can read off that the 
poles in s are located at: 

Sjn = TT^ (2.8) 

where « G No- 

Additionally, the procedure described above allows us to evaluate the coefficients of the expan- 
sion in z. numerically which we used to again test the initial conditions of the differential equations. 



3. Results 



In order to get accurate results we keep terms up to z'". Our results agree with the previous 



numerical calculation [11] within less than 1% difference. To demonstrate the convergence of the 



power expansions, we show in Figure g the form factors defined in ( |1.3[ ) as functions of s, where 
we include all orders up to z.^, z^ and z'". We use as default value z = 0.1 such that the cc- threshold 
is located ats = 0.4. One sees from the figures that far away from the cc-threshold, i.e. for s > 0.6, 
the expansions for all form factors are well behaved. 

The impact of our results on the perturbative part of the high ^^-spectrum 

is shown in Figure ^ (left), where we used the same parameters as in [Q]. The finite bremsstrahlung 
corrections calculated in [Q] are neglected. From Figure ^ (left) we conclude that for jj. = nih the 
contributions of our results lead to corrections of the order 10% — 15%. Integrating R{s) over the 
high s region, we define 

Rhigh = C dsR{s). (3.2) 
J().6 

Figure ^ (right) shows the dependence of the perturbative part of /?high on the renormalization scale. 
We obtain 

/?high,pert = (0.43 ± 0.01 (Ai)) X 10-^ (3.3) 

where we determined the error by varying /i between 2 GeV and 10 GeV. The corrections due to 
our results lead to a decrease of the scale dependence to 2%. 
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Figure 2: Real and imaginary parts of the form factors as functions of s. To demonstrate the conver- 
gence of the expansion in z we included all orders up to z*, and z^" in the dotted, dashed and soUd lines 
respectively. We put n = mi, and used the default value z = 0.1. 
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Figure 3: Perturbative part of R{s) (left) and /?high (right) at NNLL. The solid lines represents the NNLL 
result, whereas in the dotted lines the order corrections to the matrix elements associated with Oi,2 are 
switched off. In the left figure we use jx = mi,. See text for details. 
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